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® 要 : 采用 基于 Doi 理论 的 微 -宏观 双 尺 度 模型 研究 了 液晶 聚合 物 在 均匀 剪 切 流 场 中 的 微观 结构 。 为 保 
证 模拟 结果 的 可 靠 性 ， 在 模型 求解 中 率先 使 用 了 高 精度 、 高 稳定 性 的 无 网 格 方法 ， 并 通过 对 均 
匀 流 场 中 均 质 液晶 聚合 物 系 统 的 数值 模拟 ， 验 证 了 该 方法 的 有 效 性 。 文 中 着 重 研究 了 De 数 对 平 
板 Couette 流 中 液晶 聚合 物 微观 结构 的 影响 。 模 拟 得 到 了 液晶 聚合 物 分 子 四 种 典型 的 微观 运动 结 
构 ， 即 弹性 稳 态 、 周 期 翻转 、 周 期 摆动 和 随 流 取向 。 定 量 分 析 了 De 数 对 两 种 周期 运动 的 影响 ， 
预测 了 各 种 结构 下 分 子 的 有 序 程度 和 可 能 出 现 的 缺陷 ， 并 通过 对 取向 角 的 详细 分 析 ， 发 现 了 两 种 


边界 层 效应 。 
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1 引言 


液晶 是 有 别 于 液体 和 晶体 的 一 种 独立 物质 形态 ， 它 可 以 像 液体 一 样 流动 ， 同 时 又 具有 分 子 
取向 的 某 种 有 序 性 。 正 是 由 于 这 种 有 序 性 的 存在 ， 使 得 液晶 具有 高 强度 、 高 模 量 、 低 热 涨 系数 
和 良好 的 光电 特性 等 优点 。 液 晶 的 加 工 过 程 多 包含 剪 切 流动 和 拉 伸 流动 。 研 究 液 晶 分 子 在 流 场 
中 的 运动 规律 ， 对 从 宏观 尺度 控制 和 保持 液晶 分 子 的 高 度 有 序 排列 具有 重要 的 现实 意义 。 

根据 生产 条 件 的 不 同 ， 液 晶 被 分 为 热 致 液晶 和 溶 致 液晶 。 通 过 加 热 形成 液晶 相 的 称 为 热 至 
液晶 ， 通 过 加 入 溶剂 形成 液晶 相 的 称 为 溶 致 液晶 。 目 前 的 研究 主要 集中 在 溶 致 液晶 。 近 几 十 
年 来 ， 不 少 学 者 已 在 这 方面 做 出 了 一 些 突出 的 贡献 。 首 先 对 刚性 长 棒状 高 分 子 液 晶 溶 液 做 出 
理论 解释 的 是 美国 物理 学 家 Onsager!!! 和 美国 高 分 子 学 家 Floryl 引 。 在 小 分 子 液晶 方面 ， 有 著名 
的 Maier-Saupe 平 均 场 理论 引 。 在 数值 模拟 方面 ， 学 者 使 用 Doi 模 型 身 (适用 于 稳 态 流 场 ) 成 功 
模拟 出 液晶 聚合 物 溶液 众多 复杂 的 分 子 动力 学 现象 。 

平板 Couette 流 是 最 简单 也 是 最 常见 的 剪 切 流动 ， 整 个 流 场 拥 有 唯一 的 前 切 率 ， 属 于 均匀 
剪 切 流动 。 但 是 宏观 上 的 均匀 ， 并 不 能 简单 地 推演 到 微观 。 由 于 微观 边界 效应 的 存在 (边界 会 
阻碍 分 子 的 转动 )， 液 晶 聚 合 物 的 微观 结构 在 整个 流 场 中 并 不 均 勾 。 边 界 效应 破坏 了 分 子 取 向 
的 有 序 性 ， 从 而 会 形成 缺陷 。 

本 文采 用 基于 Doi 理 论 的 微 -宏观 双 尺 度 模 型 来 研究 平板 Couette 流 中 液晶 聚合 物 浓 溶 
液 的 微观 结构 。 模 拟 中 需要 直接 求解 描述 分 子 取向 扩散 的 Smoluchowski 方 程 由 ， 求 解 的 
精度 和 稳定 性 对 模拟 结果 具有 重要 影响 。 因 而 模拟 过 程 对 数值 方法 的 精度 和 稳定 性 具有 
较 高 要 求 。 无 网 格 方法 56] 是 近年 来 出 现 的 一 种 新 型 数值 方法 ， 无 单元 Galerkin (Element 
Free Galerkin, EFG) 方法 中 是 无 网 格 方法 的 一 种 。 该 方法 在 构造 近似 函数 时 ， 只 需要 节点 
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的 位 置信 息 ， 并 且 计 算 精 度 高 、 稳 定性 好 。 为 了 保证 模拟 结果 的 可 靠 性 ， 本 文 在 真 接 求 
解 Smoluchowski 方 程 和 宏观 流 场 时 均 采 用 EFG 方法 。 

文中 通过 研究 De 数 对 模拟 结果 的 影响 ， 得 到 了 平板 Couette 流 中 液晶 聚合 物 的 四 种 典型 微 
观 结构 。 并 分 析 了 De 数 与 周期 运动 的 关系 ， 研 究 了 分 子 有 序 程度 的 变化 ， 预 测 了 可 能 出 现 的 
液晶 缺陷 ， 指 出 了 边界 层 效应 的 存在 。 


2 ”液晶 流体 模拟 的 双 尺 度 模型 


液晶 聚合 物 溶液 可 视 为 流动 过 程 满足 Navier-Stokes 方 程 的 连续 介质 。 为 了 能 够 全 面 揭示 液 
晶 分 子 在 流 场 中 的 微观 结构 ， 本 文 的 数学 模型 在 大 分 子 尺 度 上 考虑 了 流 场 对 液晶 分 子 取向 扩散 
运动 的 影响 ， 同 时 在 宏观 尺度 上 考虑 了 液晶 分 子 在 流动 过 程 中 对 流 场 产生 的 应 力 。 

2.1 液晶 分 子 取向 扩散 方程 

假设 液晶 分 子 具有 棒状 结构 ， 设 棒 的 长 度 1 远 大 于 直径 6b，v 为 棒状 分 子 的 数目 密度 。 采 用 
概率 密度 函数 f (x, u, t) 来 表示 t+ 时 刻 质 心 位 于 x 的 硬 棒 指 向 平行 于 方向 的 可 能 性 ， 这 里 为 
单位 向 量 。 当 同时 考虑 转动 扩散 和 平 动 扩散 时 ，f(x, u, t) 将 满足 如 下 的 Smoluchowski 方 程 


+ (vf) 二 FY {Disa + Did— uw) - (Yf + /70)} 


十 


Dr 2. (Rf + fRU)—R-(uxn-uf), (1) 
kgT 

这 里 KB 为 Boltzmann 常数 ， 人 为 热力 学 温度 ，v 为 宏观 流 场 速 度 ，k = (Vv)? 为 速度 梯度 张 
R, R=ux 6/6u 为 转动 算 子 ， 了 ,为 旋转 扩散 系数 ， 刀 | ADL 分 别 为 平行 和 科 直 于 取向 方向 
的 平 动 扩 散 系数 ， 并 且 近 似 满 足 万 | = 2D1，U 为 排斥 体积 势 ， 常 取 为 Maier-Saupe 势 8] 


U= UokeT | ju x ul f(u')du’, (2) 
Iu’|=1 
Lo 是 与 zb12 成 正比 的 无 量 纲 参数 ， 用 以 反映 聚合 物 溶液 的 浓度 。Uo = 0 表示 稀 溶 液 ，Uo > 
0 表示 浓 溶液 。 液 晶 聚 合 物 应 力 可 以 通过 下 式 计算 


72, = 3vkBT Sag — v ((u x RU)aug) + Sern (UaUgURUL) (3) 
其 中 a, 8, k, l= 1 2,3，(.) 表 示 对 于 分 子 指向 的 系 综 平均 ， 为 分 子 旋转 摩擦 系数 ，5a8 = 
(uaug 一 6a8/3)，6a8 为 Kronecker 符 号 。 

2.2 ”宏观 流 场 控制 方程 

假设 液晶 聚合 物 溶液 具有 不 可 压缩 的 性 质 ， 则 流 场 控制 方程 为 四 


V-v=0, (4) 
vit(v-V)v =—-Vp+nsAvt V-7? 4+ F*, (5) 


v 为 流 场 速度 ，p 为 运动 学 压力 ，7s 为 溶剂 的 动力 学 粘度 ，7z2 为 液晶 聚合 物产 生 的 粘 弹性 应 
Jj, Fe 为 由 于 分 子 长 程 作用 而 产生 的 体力 。 
2.3” 微 -宏观 双 尺 度 模型 的 无 量 纲 化 


设 Lo 为 流动 区 域 特征 长 度 ，W 为 流动 特征 速度 ，To = Lo/Vo 为 流动 时 间 ， 此 时 可 将 微 - 宏 
观 烟 全 方程 无 量 纲 化 为 


V-v=0, (6) 
vit(v-V)v = -Vp+ Save gv, T? +F°), (7) 

Hf 9p) = Ev. {D'a un) (Vf + V0) 
+ FR (Rf + fRU) -R (ux ruf), (8) 
U= vo f ju x u’|? f(u’)du’, (9) 

lu’l=1 
TBa = 3Sap — ((u x RU)aug) + nn (uauguku), (10) 
(VU), (11) 
其 中 
_ Vo/Lo l _ Volo o Der Boa’ In(i/b) 
人 


3 ”无 网 格 方法 及 其 有 效 性 验证 


本 文 运 用 EFG 方法 求解 液晶 流体 模拟 的 双 尺 度 模 型 。 该 方法 采用 移动 最 小 二 乘 (Moving 
Least Squares, MLS) 逼近 四 近似 函数 ， 不 需要 对 区 域 进行 网 格 剖 分 。 
3.1 移动 最 小 二 乘 逼近 
在 区 域 Q 中 ， 未 知 标量 函数 v(x) 的 移动 最 小 二 乘 逼近 可 以 表示 为 
= 》 p;(x)a;(x) = p’ (x)a(x), (12) 
j=l 
ep aA AIA IEE, pl (x) 为 基 函 数 向 量 ，mm 为 基 函 数 个 数 。 为 了 确定 a(x)， 定 义 上 
述 逼 近 误 差 的 加 权 Le 范 数 为 


N 
三 Zua — x;)[p7(x:)a(x) — u], (13) 
式 中 N 为 节点 总 数 ，u = u(x) 为 函数 在 节点 xi 处 的 函数 值 ，w(x 一 xi) 为 节点 权 函 数 。 为 保 


iE u(x) 为 最 佳 近似 ， 需 使 R(x) 达到 最 小 ， 即 满足 
OR 
oat A(x)a(x) — B(x)U = 0, (14) 
这 里 
N 
= (unuz: un)”, A(x) = > w(x — x<)p(x:)p7 (xi), 


B(x) = (w(x ~ x1)P(x1), w(x — x2)p(X2),.… ,w(x — xn )P(xw)), 
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从 (14) 式 中 求 得 a(x) 并 代入 (12) 式 中 ， 得 到 函数 w(x) 的 MLS 逼近 为 
w(x) = p7 (x)A7 t BU = T(x)U, (15) 


JEP BT (x) = [y1(z), y2(z),--+ ,ON(zZ)]。 

3.2 有效 性 验证 

Smoluchowski 方 程 的 求解 精度 对 液晶 聚合 物 的 模拟 结果 具有 重要 影响 。 当 使 用 有 限 差分 
法 (Finite Difference Method, FDM) 计算 时 ， 误 差 累计 将 会 破坏 概 举 密 度 函 数 的 守恒 性 ( 归 一 
ME) 或 者 出 现 负 的 取向 概率 密度 值 ， 即 


27 


f(9,t)d9 关 1， 或 3f(0,t) <0, (16) 


为 了 说 明 EFG 方法 可 以 很 好 地 避免 上 述 问 题 ， 下 面 给 出 均匀 流 场 中 均 质 系统 的 模拟 结果 。 
3.2.1 均匀 流 场 中 均 质 系统 模拟 的 数学 模型 
对 于 均 质 系统 ， 取 向 概率 密度 函数 与 空间 位 置 无 关 ， 因 而 在 Smoluchowski 方 程 中 只 需 考 虑 
转动 扩散 。 此 时 (8) 式 简化 为 


oF (vf) = ER- (Rf +fRU)-R: (ux ruf), (17) 

f 与 空间 位 置 无 关 ， 且 流体 不 可 压 ， 因 而 了 :vf) =V:Yf+fY:v 三 0，(8) 式 进一步 简化 为 
ð 

Of LLR. (Rf +fRU)-R-(uxr-uf), (18) 


所 以 对 均匀 流 场 中 均 质 系统 的 数值 模拟 ， 只 需要 在 流 场 一 点 处 求解 方程 (18)。 当 分 子 的 转动 限 
定 在 平面 内 ， 并 且 考 虑 Maier-Saupe 势 时 ， 方 程 (18) 具体 为 


0 ca i Z f 
有 = 0% (f+ fUo Í sin 2(0 — 6") f (0, t)d6" ) 


—Og [ (kalcos20 + k22 cos@ sinb — «k11 cos@sin — Ky2sin76) f ] : (19) 


3.2.2 EFG 离散 格式 
当时 间 方 向 采用 Crank-Nicolson 格 式 离散 时 ， 在 第 n 十 1(n = 0,1,2,3,.…) 个 时 间 步 ， D 
程 (19) 的 EFG 离散 格式 为 


N N 
元 (2er) 本 saz (Oris X does fh 


j=l 


Qn aul 
+ 而 人 sin2(6— 0°) (6,0")d0" > 57") 
0 


了 一 1 
i N 
= 5 ( ov. | 210080 + (K22 — K11) cos @ sin 6 一 K128in"6 | 2 wfr) 
三 5{ > 二 (8 vi, fe + f°U ETG £")d6” ) 
2 De OPi: Jp 0 0 


+( doi; [k21 cos? 6 + (Kk22 — K11) cos 0 sin ô — k12 sin? 0) f” ) i: (20) 
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其 中 i = 1,2,:… ,N。 求 解 时 ， 给 定 周期 边界 条 件 ， 并 采用 Lagrange 乘 子 法 来 施加 边界 条 件 。 
3.2.3 ”方法 检验 
24 


时 ， 方 程 (19) 有 精确 解 


ft(0) = Fp ( oO cos? )， (21) 


常数 C 可 以 通过 取向 概率 密度 函数 的 守恒 性 来 确定 。 当 De = 1 时 ，C = 15.238824. € 
义工 范 数 误差 、Za 范 数 误差 和 失 恒 率 ” 分 别 为 


t h 25 t h 
sree Fle, = f 1O -O lae, 
Qn / 
| Ro fll, = (| | ft(0) “J Poan) 


"=|/ “pa -1|/ f “eeeag=| f " hoaa- 1, 


其 中 ft ARER. De = 1 时 ， 分 别 采 用 EFG 方法 和 Crank-Nicolson 格式 的 FDM 对 模型 进行 
求解 。 计 算 结果 的 比较 如 表 1 和 表 2 所 示 。 从 表 1 和 表 2 可 以 看 出 ， 与 FDM 相 比 ，EFG 方法 
具有 较 高 的 计算 精度 ， 可 以 很 好 地 保证 取向 概率 密度 函数 的 守恒 性 ， 从 而 说 明了 EFG 方法 的 
有 效 性 。 


表 1: 计算 精度 比较 


节点 数 


0.020735 
0.010851 
0.007360 
0.005569 
0.004481 


6.83915e-5 0.012508 3.70626e-5 
4.33777e-5 0.010064 2.32182e-5 
RQ: 失 恒 率 比较 


í Soe | 41 81 121 
EFG 2.611e-16 1.295e-15 1.173e-15 
0.01859 0.00995 0.00678 


81 


3.237e-15 
0.00514 


第 2 期 王晓东 等 ;均匀 剪 切 流 场 中 液晶 聚合 物 微观 结构 的 无 网 格 模拟 199 
a I I ET ORE AS AS PE Ue A GORI oe 


4 平板 Couette 流 中 双 尺 度 模 型 的 无 网 格 离散 


4.1 平板 Couette 流动 
对 于 平板 Couette 流动 ，z 轴 的 任意 垂直 截面 上 流动 情况 相同 ， 即 速度 和 压力 不 随 z 的 变化 
而 变化 。 此 时 方程 (6)-(11) 具体 为 


T= 
Me = “Rent + A (22) 
Of e? / 2 
Be = pg» PLC + sin?O)(fy + fUy)] + B60(fo + fUo) + Oo(fuvsin? 8), (23) 
U(y,8) = Uo ie sin?(0 — 6) f(y, 0", t)40", (24) 
0 
D 3 De 
Th = (sin 20) — (Ugsin?0) + “ty (sin? bcos20)， (25) 


其 中 0 为 u 和 z 轴 之 间 的 夹 角 。 计 算 时 ， 先 通过 方程 (22) 求解 流 场 ， 再 通过 (23) 到 (25) 求解 应 
力 ， 不 断 重 复 ， 直 到 终止 条 件 满足 。 

实际 中 ， 区 域 边界 对 分 子 的 转动 具有 很 大 影响 ， 使 分 子 以 一 定 的 方向 “ 锚 定 ”在 边界 上 ， 
因而 如 应 该 包含 边界 效应 由 。 当 分 子 的 错 定 方向 与 x 轴 正 向 的 夹 角 为 p 时， 边界 处 取 值 修正 为 


2r 
Uw =U» f sin? (8 — 6) f(y, 0’, t)dO’ + Uosin?(8 — p), y=0,1, (26) 
0 
Ce ee ne eee 引入 有 序 参量 s(0 < s < 1) 和 取向 角 约 (一 /2 < 


$ 7/2), BUM! 
=vVa?+b?, ¢=arcsin (b/ 2s? + 2as), 


JE fa = (cos26)/(1), b= (sin20)/(1), sHKHROT MH AFRE, s = 0 对 应 完全 无 
序 ，s = 1 对 应 所 有 分 子 一 致 取向 ; $ 表 示 取 向 方向 与 流动 方向 的 夹 角 。 

4.2 EFG 离散 格式 

由 于 分 子 的 取向 概率 密度 函数 既是 取向 空间 上 的 函数 又 是 物理 空间 上 的 函数 ， 所 以 直接 求 
解 (23) 式 给 出 的 Smoluchowski 方 程 相对 困难 。 为 了 降低 求解 难度 ， 本 文采 用 Lozomskil19 提 
出 的 时 间 分 离 技术 求解 Smoluchowski 方 程 。 该 方法 在 第 mn + 1(n = 0,1,2,3,…) 个 时 间 步 ， 
将 (23) 式 分 解 为 如 下 两 个 方程 


LGF- F 地 (fo + fUo) + ol fuysin?ð), 0 € [0,27], 20 
payee. D in26 U, 0.1 28 
Ae TF „Í 1 (1 +sin*6)(fy + fUy)], y € [0,1], (28) 


计算 时 ， 先 在 物理 空间 点 上 求解 方程 (27)， 然 后 再 在 取向 空间 点 上 求解 方程 (28)。 
记 物 理 空间 上 的 MLS 形 函数 为 pi(y) (i = 1,2,… ,Ni)， 取 向 空间 上 的 形 函 数 为 女 (9) (7 = 
+, No). EP N, No 分 别 为 物理 空间 和 取向 空间 的 离散 点 数 。 定 义 (,.)y 和 (Jo 分 别 为 
物理 空间 和 取向 空间 上 的 函数 内 积 ， 即 


1 2r 
(u,v)y = f u(yo(u)dy, (0h) = i a(9)h(6)a0 
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当时 间 方 向 采用 Crank-Nicolson 格式 离散 时 ， 双 尺度 模型 的 EFG 全 离散 格式 如 下 


(Wi, f)e + ape fo + fUo)o + 全 (oo fuysin?6)¢ 


= (pi f” Jo a (dois fo + f"Ue)o a AT Bopi, fuysin?O)o, i= 1,2,- ‚N2, (29) 
n+1 Ate? 20)( Fn+1 n+l 
(pif uy + Se Spe (Ayes, D'a + sinh) (fyi + SU), 
ft _ zI zr ， 
= (p;,f)y— Spe (PD 1(1+ sin? 0)(fy + f'Uy)),, j=1,2, ,NM (30) 
(pr u” t!)y + At 二 一 了 (8 prs uyt )y 


Re 


1] 一 
= (pr U”)y o5 e ~ At Tepe OH Ti k= 1,2, as Mi, (31) 


其 中 函数 
F = f(y;,0), f” = f(y;,9,t”), f = fly, 4), 
FH = f(y, i,t"), ut = u(y, tt), u” = uly, t”) 


及 其 对 应 的 导数 for HR fy fet, ug, dl a mt 
取向 概率 密度 函数 关于 取向 角 具 有 周期 性 ， 因 而 求解 方程 (27) 时 需要 给 定 周期 边界 条 件 。 另 
外 ， 为 了 保证 取向 概率 密度 函数 的 守恒 性 ， 在 求解 方程 (28) 时 采用 如 下 的 守恒 边界 条 件 四 


5 ”模拟 结果 


模拟 过 程 中 取 锚 定 方向 为 z 轴 方向 ， 即 p = 0; 平 动 扩散 系数 Di = 0.1; 分 子 相 对 长 度 = 
0.02; 排斥 体积 势 强度 系数 mo = 10; 上 板 驱 动 速度 为 1， 粘 度 比 7 = 0.1; Reynolds 数 Re = 
1。De 数 是 描述 聚合 物 分 子 松弛 时 间 刁 宏观 特征 时 间 之 比 的 无 量 纲 参数 ， 在 2.3 节 中 将 其 定义 
为 分 子 旋转 扩散 时 间 与 宏观 特征 时 间 的 比 。 下 面 分 析 De 数 对 均匀 剪 切 流 场 中 液晶 聚合 物 微观 
结构 的 影响 。 

5.1 四 种 典型 微观 结构 

5.1.1 弹性 稳 态 (Elastic-stability) 结构 

当 De 数 充分 小 时 ， 液 晶 聚 合 物 在 长 程 的 弹性 作用 下 保持 稳 态 ， 称 为 弹性 稳 态 结构 ， 见 
图 1。 这 种 结构 下 ， 分 子 为 了 平 衡 流 场所 引致 的 粘性 抵 ， 取 向 点 j 流 场 方向 存在 一 个 很 小 的 负 向 
偏 角 。 对 比 图 1(a) 和 图 1(b) 的 结 a 随 着 De 数 的 增 大 ， 该 偏 角 也 随 之 增 大 。 

5.1.2 周期 翻转 (Tumbling) 结 

当 De 数 逐 渐 增 大 时 ， es eRe eee. 在 非 零 转动 矩 的 作用 下 ， 液 晶 
聚合 物 分 子 开 始 出 现 周期 性 的 翻转 运动 。 图 2 (a) 是 弹性 稳 态 向 周期 翻转 的 过 渡 ， 此 时 已 经 可 
以 观察 到 周期 性 翻转 运动 ， 但 流 场 内 部 各 点 处 分 子 的 翻转 并 不 完全 一 致 。 随 着 De 数 进 一 步 增 
大 ， 流 场 内 部 各 点 处 分 子 的 翻转 趋 于 一 致 ， 液 晶 聚 合 物 出 现 稳定 的 周期 翻转 结构 ， 见 图 2(b)。 

5.1.3 周期 摆动 (Wagging) 结构 
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随 着 De 数 继续 增 大 ， 分 子 的 转动 开始 受到 分 子 问 排斥 体积 作用 的 抑制 ， 运 动 形式 由 周期 番 
转向 周期 摆动 转变 ， 见 图 3。 在 这 种 结构 里 ， 分 子 绕 流 场 方向 不 断 摆动 。 并 且 随 者 De 数 的 增 
大 ， 分 子 转动 的 难度 逐渐 增 大 。 


(a) De = 0.01 (b) De=0.1 
图 1: 弹性 稳 态 结构 的 取向 角 演 化 图 


(a) De = 0.2 (b) De = 2.0 
图 2: 周期 翻转 结构 的 取向 角 演 化 图 


(a) De = 15 (b) De = 18 
图 3: 周期 摆动 结构 的 取向 角 演 化 图 


5.1.4 随 流 取向 (Flow-aligining) 结构 

X De 数 足够 大 时 ， 分 子 的 转动 完全 被 抑制 ， 整 个 流 场 再 次 早 现 稳 态 结构 ， 所 有 分 子 一 致 沿 
流 场 取向 ， 称 为 随 流 取向 结构 。 图 4(a) 给 出 了 周期 摆动 向 随 流 取向 转变 的 过 渡 状 态 ， 此 时 分 子 
依然 在 周期 摆动 ， 但 是 摆动 幅度 随 着 时 间 的 推进 不 断 衰 减 。 图 4(b) 给 出 的 则 是 分 子 运动 完全 进 
入 随 流 取向 结构 的 状态 。 


(a) De = 20 (b) De = 50 
图 4: 随 流 取向 结构 的 取向 角 演 化 图 


5.2 De 数 对 周期 运动 的 影响 

5.1 节 指出 ， 液 晶 聚 合 物 在 均匀 剪 切 流 场 中 存在 两 种 具有 周期 性 的 微观 结构 ， 即 周期 翻转 结 
构 与 周期 摆动 结构 。 图 5 给 出 了 De 数 对 这 两 种 结构 的 影响 。 从 图 5(a) 可 以 看 出 ， 在 进入 周期 
翻转 运动 之 初 ， 翻 转 周 期 随 De 数 的 增 大 ， 明 显 减 小 ， 随后， 翻转 周期 达到 下 限 ，De 数 对 其 
不 再 具有 影响 。 从 图 5(b) 可 以 看 出 ， 在 整个 周期 摆动 阶段 ， 摆 动 周期 随 De 数 的 增 大 而 减 小 。 
与 翻转 阶段 相 比 ， 摆 动 阶段 周期 的 变化 比较 缓慢 。 图 5(c) 给 出 了 分 子 周 期 摆动 时 ， 摆 动 幅度 
Bi De 数 的 变化 。 从 图 中 可 以 看 出 ， 正 向 ( 逆 时 针 方向 ) 摆动 幅度 始终 小 于 反 向 ( 顺 时 针 方向 ) FE 
动 幅度 。 两 者 均 随 De 数 的 增 大 而 减 小 ， 并 且 差 中 逐渐 增 大 。 


(RE) 


014 15 16 17 18 
De 


(a) 翻转 周期 (b) 摆动 周期 (c) 摆动 幅度 


图 5: De 数 对 周期 运动 的 影响 
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5.3 ”有 序 程度 及 缺陷 预测 

研究 分 子 取向 的 有 序 性 并 预测 可 能 存在 的 缺陷 ， 是 液晶 聚合 物 数值 模拟 的 重要 内 容 。 
图 6 给 出 了 几 组 有 序 参 量 随时 间 的 演化 等 值 云图 。 图 6(a) 对 应 弹性 稳 态 结构 。 从 图 6(a) 中 可 以 
看 出 ， 有 序 参 量 取 值 保持 在 0.8 以 上 ， 边 界 附近 存在 一 些 带 状 缺 陷 。 图 6(b) 对 应 弹性 稳 态 结构 
向 周期 翻转 结构 的 过 渡 阶 段 。 此 时 ， 有 序 参量 在 流 场 大 部 分 区 域 保持 在 0.8 以 上 ， 边 界 附 近 除 
了 存在 带 状 缺陷 ， 还 出 现 了 明显 的 点 状 缺 陷 。 图 6(c) 和 图 6(d) 对 应 周期 翻转 结构 。 有 序 参量 
出 现 了 较 大 幅度 的 周期 性 变化 ， 取 值 周期 性 地 浮动 于 0.5 到 1.0 之 间 ， 边 界 附近 的 带 状 缺 陷 逐 
渐 消失 ， 点 状 缺陷 依然 比较 显著 。 图 6(e) 和 图 6(f) 对 应 周期 摆动 结构 。 有 序 参量 变化 的 幅度 
较 前 有 所 增加 ， 取 值 在 0.1 到 0.9 之 间 周 期 性 变化 ， 边 界 附近 的 点 状 缺 陷 程度 较 前 有 所 减弱 。 
图 6(g) 和 图 6(h) 分 别 对 应 周期 摆动 向 随 流 取向 的 过 渡 阶 段 和 随 流 取 向 阶段 。 此 时 ， 有 序 参量 
不 再 周期 性 浮动 ， 取 值 保持 在 0.7 到 0.8 之 间 ， 点 状 缺 陷 消失 ， 在 液晶 聚合 物流 体 完全 处 于 随 流 
取向 状态 时 ， 模 拟 结果 无 法 预测 任何 缺陷 。 

a 

01 02 03 04 05 06 07 08 0.9 1.0 


i 0.1 0.2 03 04 05 06 07 0. 


0 7 190 ~ 150 "200 0 50 100 T50 200 


1 0.1 0.2 03 04 05 06 0.7 08 09 1.0 


(g) De = 30.0 (h) De = 50.0 


图 6: 有 序 参 量 随时 间 的 演化 等 值 云图 
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5.4 De 数 与 边界 层 效 应 的 关系 i 

由 于 边界 会 阻碍 分 子 的 转动 ， 因 而 在 周期 翻转 结构 或 周期 摆动 结构 中 可 以 观测 到 明显 的 边 
界 效 应 。 图 7 给 出 了 几 组 取向 角 随 时 间 的 演化 等 值 云图 ， 其 中 等 值 线 不 同 程度 的 弯曲 说 明了 流 
场 中 存在 边界 层 效 应 ， 并 且 这 种 效应 与 De 数 的 大 小 有 关 。 

为 了 清楚 地 观测 这 种 边界 层 效 应 ， 并 研究 其 随 De 数 的 变化 规律 ， 图 8 给 出 了 一 个 周期 内 取 
向 角 随 时 间 演 化 的 局 部 放大 等 值 云图 。 从 图 8(a) ~ 图 8(c) 可 以 看 到 ， 发 生 周期 运动 时 边界 层 附 
近 分 子 较 流 场 内 部 分 子 有 一 段 时间 的 延迟 ， 我 们 将 其 称 作 边 界 汪 后 效应 。 并 且 随 着 De 数 的 增 
大 ， 边 界 滞 后 效应 具有 明显 的 减弱 趋势 。 当 De 数 增 大 到 一 定 程度 之 后 ， 边 界 灌 后 效应 消失 ， 
随 之 出 现 的 是 边界 前 驱 效应 ( 见 图 8(d) 必 图 8(f))。 出 现 边界 前 驱 效应 时 ， 千 期 运动 中 边界 层 附 
近 分 子 较 流 场 内 部 分 子 提前 发 生 转动 。 并 且 随 着 De 数 的 增 大 ， 边 界 前 驱 效应 显著 增强 。 


(e) De = 12.0 (f) De = 15.0 


图 7: 取向 角 随 时 间 的 演化 等 值 云图 
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55 56 57 59 59 60 61 62 63 


(b) De = 2.0 


51 52 53 54 5 56 57 


(c) De = 5.0 (d) De = 10.0 


51 52 53 54 pi 56 57 58 51 52 53 54 5 56 57 58 59 60 


(e) De = 12.0 (f) De = 15.0 


图 8: 取向 角 随 时 间 演 化 局 部 放大 等 值 云图 


6 结论 


首次 运用 EFG 方法 研究 了 均匀 剪 切 流 场 中 液晶 聚合 物 的 微观 结构 。 该 方法 很 好 地 保证 了 取 
向 概率 密度 函数 的 计算 精度 及 守恒 性 ， 从 而 确保 了 模拟 结果 的 可 靠 性 。 模 拟 研究 表明 : 

1) 随 着 De 数 的 增 大 ， 液 晶 聚 合 物 会 依次 呈现 出 四 种 典型 的 微观 结构 ， 即 弹性 稳 态 结构 、 
周期 翻转 结构 、 周 期 摆动 结构 和 随 流 取 向 结构 。 

2) 在 周期 翻转 结构 中 ， 随 着 De 数 的 增 大 ， 翻 转 周 期 起 初 迅速 减 小 ， 后 来 趋 于 一 稳定 值 ; 
而 在 周期 摆动 结构 中 ， 随 着 De 数 的 增 大 ， 摆 动 周期 不 断 缓慢 减 小 ， 并 且 摆 动 幅度 不 断 衰减 。 

3) 四 种 结构 中 ， 分 子 都 能 保持 较 高 的 有 序 度 。 并 且 受 分 子 周 期 运动 的 影响 ， 有 序 度 也 会 
呈现 周期 性 变化 。 

4) 除了 随 流 取向 结构 外 ， 在 其 他 三 种 结构 中 不 同 程度 地 预测 到 了 液晶 缺陷 。 其 中 ， 在 弹 
性 稳 态 结构 中 会 出 现 带 状 缺陷 ， 周 期 翻转 结构 中 会 出 现 明 显 的 点 状 缺 陷 和 轻微 的 带 状 缺 陷 ， 周 
期 摆动 结构 中 只 有 轻微 的 点 状 缺 陷 。 

5) 液晶 聚合 物 处 在 周期 翻转 结构 或 周期 摆动 结构 时 ， 边 界 层 附近 可 以 观测 到 严重 的 边界 
层 效 应 。De 数 较 小 时 ， 表 现 为 边界 滞后 效应 ， 并 且 效 应 随 着 De 数 的 增 大 而 减弱 ，De 数 较 大 
时 ， 表 现 为 边界 前 驱 效应 ， 并 且 效 应 随 着 De 数 的 增 大 而 增强 。 
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Simulation of Microstructure for Liquid-crystalline Polymers in 
Homogenous Shear Flow by Meshfree Methods 


WANG Xiao-dong, OUYANG Jie, SU Jin 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi’an 710129) 


Abstract: The micro-macro dual-scale model based on the Doi theory has been used to study the 
microstructure of liquid-crystalline polymers in homogenous shear flow. The element free Galerkin 
(EFG) method with high accuracy and good stability is used to solve the model. The validity of the 
EFG method has been verified through the simulation of homogenous liquid-crystalline polymers in 
uniform flow. Four typical microstructures of liquid-crystalline polymers in plane couette flow, say 
elastic-stability, periodic tumbling, periodic wagging and flow-aligning, have been found by observing 
the affect of the De number. Moreover, we analyze the influence of the De number on two periodic 
motions, predict the degree of molecular orders and implied defects, and find the two kinds of boundary 
effects. 

Keywords: shear flow; liquid-crystalline; microstructure; dual-scale; meshfree 
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